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Abstract— A numerical solution on forced convection of Alz03-water nanofluid for different volume fractions 
is investigated for laminar flow through a parallel plate with flush mounted discrete heat sources. The model 
used for nanofluid mixture is a single-phase approach and fluid properties are considered constant with 
temperature. The finite difference method is used for solutions and four different volume fractions are 
considered varying from 0% to 4%. A fully developed laminar velocity profile is considered and the parallel 
plate is assumed as heated with three discrete heat sources flush mounted to the top and bottom plate with 
the same lengths. Uniform wall temperature boundary condition is taken for discrete heaters. Peclet numbers 
are in the range of 20-100. For comparison and validity of the solution the results for a classical problem, 
laminar flow through a parallel plate which is heated at the downstream region with constant temperature, 
are obtained. Results are presented in terms of bulk temperature, heat flux, and local Nusselt number. Heat 
transfer is enhanced with the particle volume concentration. For comparison, pure water results are also 
shown in the figures. At the locations where heat is applied the heat flux values decrease as the volume 
fraction increase and the bulk temperature values are higher for the higher volume fractions at the heated 
locations. As the volume fraction increases the local Nusselt number can increase up to 30% than to pure 
water. 


Keywords— forced convection, heat transfer enhancement, nanofluid, flush mounted discrete heating, 
parallel plate, numerical investigation. 


I. INTRODUCTION There are papers concerning mixed convection heat transfer 


The classical problem of the laminar flow of fluids in the from discrete heat sources [8], [9], [10], [11], [12], [13] and 


parallel plates has been investigated by many researchers [3]. 


for various boundary conditions. A detailed survey of these 
works can be found in [1]. There are also works about 
discrete heat sources for pipe and channel flows in 
literature. A comprehensive review of work on discrete heat 
sources in plates were presented by [2] and [3]. 


Two-dimensional conjugate heat transfer from flush 
mounted heat sources in horizontal plates numerically by 
[4], [5] and [6]. experimental and numerical analysis is 
performed by [7] to determine convection heat transfer from 
a single and inline four-row array of heat sources. 
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An innovative way to enhance the heat transfer rate in 
thermal equipment is using solid nanoparticles to improve 
the fluid thermal conductivity [14]. Nanofluids constitute 
base fluids and metallic or non metallic particles having size 
less than 100nm and are known to have better thermal 
conductivity than conventional fluids, [15]. A detailed 
review of convective heat transfer enhancement can be 
found in [16] and [17]. 


There is not much paper concerning nanofluid and discrete 
heat sources in the literature. According to our 
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acknowledgement survey, [18], [19] and [20] are the papers 
that can be found in the literature. 


In this paper, forced convection of nanofluids (Al203-water) 
through two-dimensional horizontal parallel plate with 
discrete heat sources is investigated numerically. It is 
assumed as six heat sources fixed on top and bottom plate 
with the same length (Ln*=24) and with the same space to 
each other, three at bottom and three at top. The fluid 
assumed as enters to system with constant temperature, To. 
The discrete heaters have constant temperature of 7). The 
velocity profile is considered fully developed far away from 
the heating region. Three different values of Peclet 
numbers, 20 50 and 100 are taken and Al20O3-water is 
considered as nanofluid. Four different volume fractions as 
1%, 2% 3% and 4% are taken for solutions and for 
comparison 0% (pure water) results are given. The variation 
of bulk temperature, heat flux, and local Nusselt number 
values according to the axial position are presented 
graphically. 


Il. MATHEMATICAL MODEL 
2.1. Formulation of the Problem 


In the present work, a steady heat transfer problem is 
investigated in two-dimensional, horizontal parallel plates 
for thermally developing laminar flows, considering fluid 
axial conduction. The schematic diagram and the coordinate 
systems are given in Figure 1. The parallel plates are 
considered infinite in length and two regional. At far 
upstream (x*=-0o) fluid enters the channel with a uniform 
temperature of T,. At the downstream of the channel, it is 
assumed that three discrete heaters mounted both up and 
down plates symmetrically which are maintaining constant 
wall temperature, 7). The length of the heated and unheated 
locations is considered the same length. Totally six, three of 
them placed bottom and the other three heaters placed top 
plate symmetrically. The lengths of the discrete heaters are 
taken as Lr, and the space (unheated length) are equal to this 
value. Viscous dissipation is neglected at the solution of the 
energy equation. For the nanofluid it is assumed that; 


-there is a thermal equilibrium between the base fluid (pure 
water) and the nanoparticles. 


-the nanofluid flow is laminar and steady. 
-the nanofluid is Newtonien and incompressible [19]. 


- the single-phase model with physical and thermal 
properties, all assumed to be constant with temperature. 
[14]. 
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T=T, x0 T} To T, To T, 


Figure 1. Schematics of the problem and the coordinate 
system 
The fluid is water based nanofluid containing Al203 
nanoparticles. The solid volume fractions, ọ, have been 
varied between 0% to 4%. The thermophysical properties of 


the base fluid and nanoparticles are tabulated in Table 1. 


Table 1 Properties of nanofluid [19] 


Property Water Alumina 
(Al203) 
i Dkg/m3] 997.1 3970 
Cp [I/kgK] 4179 765 
k [W/mK] 0.613 40 
u[Pa s] 998e-6 - 


The flow is assumed developed from the channel entrance. 
By analytical solution for the parallel plate the velocity 
profile is given in Eq. 1. 


2 
u(y) = w|2-(2) | 
(1) 


Under these assumptions the described problem can be 
formulated in non-dimensional form as follows. The energy 
equation is 


6Pe,,(y*-y# 2 a8 (e2) 2 G - 
Oy * 


ox* ox* Ox * 
(2) 

Boundary conditions; 
at x* = —oo 0=0; (3.a) 

00 
at y*=0.5 ——=0; (3.b) 

Oy ok 
at y*=0 and y*=1 ;x*>0 0 = 1 (at wall with 
discrete heat sources) (3.c) 


at y* = 0 and y* =i 7 x*>0 0 = 0 (at wall without 


discrete heat sources) (3.d) 


Oy * 
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at y*=0 and y*=1; x*<0 O=O0 (the upstream 


region) (3.e) 
Non dimensional parameters for the problem are defined 
as; 

me iS x2 L*= h yea 

T=. d "ad d 
du, P C p k, 
Pey=——*, kta, (4) 
Kapp koy »b 


Dimensionless fluid bulk temperatures, 6,, heat flux 
values, Vy , and local Nusselt numbers, Nu, can be 


calculated as follows: 


1 
0, =6|(y*-y*)O, dy*; (5) 
0 
._ (00) 
dy ma ôy * > (6) 
Nie ie, (7) 
0, — 0, 


2.2. Thermophysical properties of nanofluids 


In the present analysis, the nanofluid flow is considered as 
a single phase incompressible flow. The simplest approach 
for the single phase assumption is the direct usage of the 
governing equations of pure fluid flow with the 
thermophysical properties of the nanofluid in consideration. 
[16]. Nanofluids may be considered as Newtonian fluids for 
low volume fractions, such as up to 10% and for small 
temperature increases. 


The density of the nanofluid based on nanofluid volume 
fraction is given by [21] as follows. 


Py =(1-9)p; +9P, 


(8) 
For a nanoparticle liquid suspension the parameter 
(pCp) nf Of the nanofluid is given as [22] 
(oc, ), =(1-ø\pc,) + olpc, ), a 


Nanofluid viscosity is one of the important thermophysical 
parameter due to its direct effect on pressure drop in forced 
convection. A detailed survey on viscosity of nanofluid can 
be found in [17] The effective dynamic viscosity of the 
nanofluid is given by [23] as 


My = L,(1+39.119+ 533.99") 


nanoparticles) (10) 


(for AlzO3 
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Many researchers suggest using nanofluids due to its higher 
thermal conductivity than to pure base fluid. A detailed 
literature review can also be found in [17]. In order to take 
the effect of thermal dispersion into account, the energy 
equation for the forced convection heat transfer of pure 
fluids is modified by replacing the thermal conductivity 
terms with an effective thermal conductivity which is 
defined as follows [18]. 


kos = Ky + Ka (11) 


Kea =Clpc,) ugd, r, 42) 
where, k,is nanofluid thermal conductivity, and ka is called 
the dispersed thermal conductivity, which is proposed to be 
calculated by using the following expression. Here, C is an 
empirical constant that should be determined by matching 
experimental data, ux is local flow velocity, d, is particle 
diameter, and r, is tube diameter. 


2.3. Numerical Procedure 


The system of equations (2) is solved by a numerical finite 
difference approach. Once the flow pattern is known, the 
heat transfer can be calculated by a numerical solution of 
the energy equation. The conductive terms are discretized 
by central difference schemes and convective term in the 
energy differential equation by an exact method defined in 
[24]. This method of discretization is a two-dimensional 
version of the “exact or exponential scheme” defined by 
[25]. The following discretization equation is obtained for 
an interior (non boundary) nodal point (i,j). 


a; Fj = tiaj stp + aiaj aj +H ja Bj ta js, 
(13.a) 
where 
6P | * a 
On; = a (ay): 43.) 


expl Pey; ~ yp J -1 


aij = 6Pdy; yj os Ay; ) 
ex} 6Pely; = | a 1 
l (13.c) 
Gi jst =5 A); (13.d) 
j+l 
k * 
ai= a (ax: ); (13.e) 
j-l 
3 


i,j-l 
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a; j = Gi j +44; +a 


i 1+) 44. (13.f) 


t, J+ 
Due to symmetry, the solution domain is bounded between 
the wall and the central line of the parallel plates. Uniform 
grid spacing is used at solutions both in x and y directions. 
The optimum number for the grid system is found to be 
160x60 in x and y directions respectively. Gauss-Siedel 
iteration technique is used for obtaining the temperature 
distribution. Convergence limit is taken to be 107. 


I. CODE VALIDATION 


A computer code has been developed to solve the energy 
equation for discrete heat sources in a parallel plate. To 
verify the reliability of the program a solution was obtained 
for a classical parallel plate problem. By changing boundary 
conditions of classical problem, the developed program is 
used for discrete heating. In this classical problem, at the 
upstream region, the inlet temperature of the fluid is 
assumed as 7,. And At the downstream uniform wall 
temperature (T;) boundary condition is taken. Also fully 
developed laminar flow is assumed for steady state 
problem. The schematic of the classical problem is given in 
Figure 2. 


u(y) = tty 
G 


x=+00 


T=T x=0 T; 


0 


Fig.2. Schematics of the problem and the coordinate 
system for code validation 


For fully developed laminar channel flow under constant 
wall boundary condition, the Nusselt number is 7.54070087 
[1]. The Nusselt numbers obtained by the computer program 
according to Peclet numbers and relative deviations from 
the value of 7.54070087 are given in Table 2. As can be seen 
from Table 2 the relative deviation between the literature 
and the obtained result is a maximum of 0.22 percent. It can 
be said that the results agree well with the literature. 


Table 2. Nusselt number values for various Peclet 


numbers. 
Peclet Nusselt Number Relative Error 
Number e (%) 
1 7.5396 0.015 
2 7.5398 0.012 
5 7.5400 0.0093 
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10 7.5399 0.0106 
20 7.5346 0.0809 
50 7.5397 0.013 
100 7.5239 0.22 


For the classical heat transfer problem to obtain the effect 
of volume fractions of nanoparticles, the axial variation of 
Nusselt numbers and bulk temperatures are calculated. The 
Nusselt number for Pe=10 is given in Figure 3. As can be 
seen from Fig.3, for small Pe number considerable amount 
of heat is transferred through the upstream region due to 
axial conduction in the fluid. 


To see the effects of volume fraction of nanoparticles, the 
values of the local Nusselt numbers at flow direction is also 
given in Fig 3 for four different particle volume fractions 
from 1% to 4% for the same boundary condition. In these 
solutions, Al203-water nanofluids are taken as the working 
fluids. Pure water shows the lowest local Nusselt number 
values. The heat transfer increases with an increase in 
nanoparticle concentration which is caused by the high 
thermal conductivity of the nanofluid. The higher local 
Nusselt values are obtained for pure water at the upstream 
region of the parallel plate. And Nusselt values decrease as 
volume fraction increases at the upstream region due to the 
higher thermal conductivity of nanofluid. This values are 
coincident with the results of [26]. 


Pe=10 
AI,O, nanofluids 


Fig.3. Axial distribution of local Nusselt number at 
classical heat transfer problem for Pe=10 


For determining the effects of volume fractions of 
nanoparticles the bulk temperatures of the classical heat 
transfer problem are given. The axial bulk temperature 
distribution for four different volume fractions 1%,2%, 3%, 
and 4% is given in Figure 4. In Fig. 4.a the results are given 
for Pe=50, and Fig.4.b is for Pe=100. Form Fig.4 it can be 
seen that the thermally developing region increases as the 
Peclet number increases due to higher axial velocity. In Fig. 
it can also conclude that the effect of axial conduction could 
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not be seen at these ranges of Peclet numbers. Despite 
handling the axial conduction term in the energy equation 
there is not any heat penetration through the upstream 
region of the fluid. The effect of nanoparticles volume 
fraction can easily be seen at higher Peclet numbers. The 
lines are more distinct from each other at Fig.b. 


9b 
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As can be seen from Figure 4. the volume fractions of the 
nanoparticle do not have a significant effect on the thermal 
entrance lengths. The bulk temperature values are higher as 
the volume fraction increases due to higher thermal 
conductivity values. 
asymptotic value of 1 as expected. 


The bulk temperature gets its 


1,04 
0,84 
0,64 Pe=100 a 
J ALO, nanofluid wes 
23 o=4% 
04-4 
0,24 
] (b) 
0,04 


Ra TA o aT Veh Ca T Ca Pee Lacan Hs 
0 10 20 30 40 50 60 70 80 90 100 
x" 


Figure 4. Axial distribution of fluid bulk temperature a)Pe=50 b)Pe=100 
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0,04 
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m 
Furthermore, to verify that the solutions are grid 


independent, results were assessed based on the generalized 
Richardson extrapolation and grid converge index (GCI), 
suggested by [27]. Sample solutions were made by using 
coarse (80x30) and fine grids(320x120) assuming the 
selected grid system as medium (160x60). A second order 
method is used in GCI analysis by taking the grid 


refinement ratio as 2.0. The quantity of interest for the 
the fluid bulk temperature @, heat 
flux, qw* and Nusselt Number, Nu, at the beginning of the 
heating section, x*=0. The parameter values used for 
calculations and the results of GCI analysis are given in 
Table 3. The computations may be assumed within the 


comparison is 


asymptotic range and no further refinement is necessary. 


Table 3. Effect of grid size on @, q*w and Nu 


axial grid 
values at down- radial fluid bulk heat flux , Nusselt 
x*=0 up stream stream grid temperature (0%) q*w Number, Nu 
ee region 

Pe=50, p=0% 
fine grid 80 280 120 0.189050 4.6496745 | 9.9644497 
medium grid | 40 120 60 0.193015 4.6542981 | 10.023364 
coarse grid 20 60 30 0.201861 4.6921733 | 10.216929 
€e between medium and coarse grids (%) 0,018945467 0,00807 0,04382 
GCI between medium and coarse grids (%) 1,036172723 0,140301 4,449492 
€ between fine and medium grids (%) 0,005877726 0,00099 0,02054 
GCI between fine and medium grids (%) 0,321466841 0,017267 2,085692 


IV. RESULTS AND DISCUSSION 


The problem handled depends on Peclet number, Pe, and 
volume fraction of nanoparticle Al203-water and For 
different volume, fractions are considered. Solutions are 
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made for Pe= 20, 50 and 100 and j =1%, 2%, 3% and 4%. 
The volume fraction values are chosen at the limit value of 
single phase solution. Results are given by heat flux values, 
fluid bulk temperature values and local Nusselt number 
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values. The dimensionless distance is Ln*=24 for whole 
solutions. Also the space (unheated lengths) is equal to 24. 
In the figures the heated axial distances are presented since 
the rest of the downstream is assumed as unheated. 


In Figure 5, variations of surface heat flux values for Pe=20, 
50 and 100 are given for four different volume fractions of 
nanofluid. In this figure, positive and negative heat flux 
values can be seen. The positive heat flux values correspond 
to heated sections of the parallel plate. The locations 
without heater of the wall show negative interfacial heat 
fluxes due to lower fluid temperature than to the wall. The 
effect of axial conduction cannot be seen due to the high 
Peclet number since there is not any heat penetration to the 
upstream region. As can be seen from Fig.5, as the Peclet 
numbers increases the heat flux values increases. 


The effect of volume fractions of nanoparticles can also be 
seen in Fig.5. The heat flux values are greater for pure water 
for heated locations. 


Volume fractions of nanoparticles do not change the values 
of heat fluxes significantly. At the locations where the heat 
is applied the heat flux values decreases as the volume 
fraction increase. But in the un-heated locations, this 
phenomenon changes vice versa. In the locations without 
the heater, the lower values of heat fluxes are for pure water 
and the values of heat fluxes are increased as the volume 
fraction increases. The value of heat flux of pure water is 
lower than the other volume fraction at unheated locations. 


And as the [| increases the heat flux values increases also. 
From Fig. 5. it can also be seen that the peak values get the 
same values this can be due to value of heating and 
unheating lengths. Especially for Pe=20 and Pe=50 (Fig 5. 
a and b) the heating and unheated lengths are sufficient for 
thermal developing but it is not for Pe=100. Because of this 
phenomenon the first peak value of heat flux is slightly 
greater than others for Fig.5.c. 


In Figure 6 distribution of bulk temperatures as axial 
distance is given. Fig. 6.a, 6.b and 6.c are given for Peclet 
numbers 20, 50 and 100 respectively. In these figures the 
heated and unheated lengths are taken as 24 and equals to 
each others. As can be seen from Figure 6, there is not any 
heat diffusion through the upstream due to axial conduction 
in the fluid. In bulk temperature curves the increment trend 
occurs during the heating location and the decrement trend 
follows it at the unheated locations. The maximum and 
minimum values are bounded between 0 and 1. The bulk 
temperature values are higher for the higher volume 
fractions at the heated locations of the parallel plate. In the 
unheated locations, the bulk temperature values get lower 
values for high volume fraction than the others. This can be 
explained due to higher thermal conductivity of nanofluids. 
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The thermal conductivity increases as the volume fraction 


increases. 
4 
34 Pe=20 
Al,O, nanofluid 
24 
14 
*F o+ 
~ increasing 
jecrensing = 
at 7 @ (4,3.2.1,0) T 
-24 
-34 
Atere 1 
O 10 20 30 40 50 60 70 80 90 100 110 
se 
5 
4] Pe=50 
ALO, nanofluid 
3] 3 
2] (b) 
14 
*F 0 
-14 
— 9=0% 
-24 9=1% 
9=2% 
3+ 9=3% 
-4-4 —— 9=4% 
544 


a RE 1 
O 10 20 30 40 50 60 70 80 90 100 110 


Pe=100 
-64 Al,O, nanofluid 


0 10 20 30 40 50 60 70 80 90 100 110 

” 
Fig.5. Variation of surface heat flux. a)Pe=20, b)Pe=50 
c)Pe=100 


As can be seen from Fig. 6.a, the bulk temperature gets its 
limit value of 1 at a short distance than the heated sections. 
This is also seen for unheated sections. In these locations 
the value gets its value of 0. This can be explained by the 
slow motion of the fluid in parallel plate. The effect of 
volume fractions of nanofluids cannot be seen at lower 
Peclet numbers distinctly. The curves of volume fractions 
get clear as the Peclet numbers increases. The smaller 
Peclet number indicates slower fluid flow in the channel. As 
the Peclet number decreases the differences between bulk 
temperature curves of nanofluid decreases due to higher 
fluid conduction. 


As can be seen from Fig. 6.c, at downstream region the peak 
values of bulk temperatures move away from boundary 
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values The top value is nearly 0.9 and the bottom bulk 
temperature value is 0.05. The rapid fluid motion causes this 
difference. The effects of volume fractions of nanoparticles 
are more significant at higher Pe numbers. 


% 


0,44 


0,24 


Pe=20 


0,0 Al,O, nanofluid 


i a i a ae a a a ee i eee ee: 
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ze 
Figure 6. Axial distrubition of fluid bulk temperature. 
a)Pe=20, b)Pe=50 c)Pe=100 


Figure 7 shows the local Nusselt number values for the 
downstream section of the parallel plate. Fig. 7.a, 7.b and 
7.c are given for Pe=20, 50 and 100 respectively. As can be 
seen from Fig. 7 at the vicinity of heating section Nusselt 
numbers get their maximum value and asymptotically 
decrease to a limit value at each un-heated locations. The 
effects of volume fractions are clearer in these three figures. 
At each Nusselt graphics, the minimum values are for the 
pure water (~=0%). the Nusselt number values increase as 
the volume fraction of nanofluid increases. 
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From Fig.7 it can also seen that the values of Nusselt 
numbers are greater for greater Peclet numbers. The local 
Nusselt number values increases as the Peclet number 
increases. This result is coherent with the results of Fig.6, 
since the differences of bulk and wall temperatures are 
higher at higher Peclet numbers. The temperature difference 
is at the denominator of Eq. 7. 


As can be seen from Fig.7.a. by using 4% of nanofluid the 
Nusselt number increases up to 25% for Pe=20. From Fig. 
7.b this percentage is up to 30% for Pe=50, and from Fig.7.c 
it is 30% for Pe=100 than to pure water. Similar Nusselt 
number increment ratios are also obtained for four different 
volume fractions for classical heat transfer problem. This 
can be seen by comparison of Fig.3 and Fig. 7. This result 
is coincident with the result of Huang (2015). The results 
show that the Nusselt numbers of Al203-water nanofluids 
are significantly higher than that of base fluid because of 
their higher thermal conductivity. 
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Figure 7. Local Nusselt numbers a)Pe=20, b)Pe=50 
c)Pe=100 
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V. CONCLUSION 


This works presents an investigation of heat transfer in a 
parallel plate heated with flush mounted discrete heat 
sources with constant surface temperature. The problem is 
handled for two regional plate, infinite in length, for which 
the upstream region is completely kept at uniform 
temperature To, which is also the inlet temperature. After 
x=0, flush mounted six heaters at temperature T; placed at 
the top and bottom plates (three for top and three for bottom) 
with same spaces while the rest of the system kept at T,. The 
solution is made numerically by a finite difference method 
for hydrodynamically developed laminar flow. Al203-water 
is assumed as flowing through the parallel plates. The 
effects of volume fractions of nanofluid were obtained for 
four different values varying from 0%-4%. The heat fluxes, 
bulk temperatures, and Nusselt numbers are obtained for 
three different Peclet numbers. A classical heat transfer 
problem, which is heated with constant wall temperature 
downstream, is solved for the validity of the computer 
program. And effects of volume fractions are also 
determined for this problem. The results obtained may be 
outlined as follows. 


> The local Nusselt number values increase as the 
volume fraction increases as expected due to higher 
thermal conductivity in the classical heat transfer 
problem. 

> The volume fractions of nanoparticles do not have a 
significant effect on thermal development length on 
the classical problem. 

> At the locations where heat is applied the heat flux 
values decreases as the volume fraction increase. And 
at unheated locations it vice-versa. 

> The bulk temperature values are higher for the higher 
volume fractions at the heated locations. 

> The Nusselt numbers of nanofluids are significantly 
higher than that of base fluid because of their higher 
thermal conductivity. It can increase the heat transfer 
by 30% than to pure water. 


Nomenclature 

a constant of discretization equation (Eq. 14) 
y normal coordinate 

d width of the channel 

u axial velocity 

Pe Peclet number 

k thermal conductivity 

L length 

T temperature 


axial coordinate 


= 
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Cp specific heat at constant pressure 
q heat flux 

T temperature 

Nu Nusselt number 

C ampirical constant 

Subscripts 

nf nanofluid 

m mean 


b bulk 

eff effective 

h heated length 
w wall 
f fluid 

p particle 

p 

a 

u 


volume fraction 


thermal diffusivity 
dynamic viscosity 
ij at nodal point i,j 
Superscripts 


*Dimensioneless quantity 


Greek symbols 

0 dimensionless temperature 
p density 

Ox axial position of difference 
oy normal position of difference 
Ax axial step size 

Ay normal step size 

€ quantity in GCI analysis 
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